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1.  INTRODUCTION 


This  report  describes  the  work  performed  under  USAF  contract  No.  F19628-7S- 
C-0167  entitled  "Investigation  of  Silicon  CCD  Delay  Lines."  The  primary  objectives  of 
the  contract  include  the  analysis,  design,  and  development  of  advanced  CCDs  with  input 
and  output  circuitry  as  well  as  delay  stages  capable  of  operation  at  100  MHz  and  above 
clock  frequency. 

A  high  speed  256  stage  CCD  delay  line  with  high  frequency  input  and  output  cir¬ 
cuits  was  designed. 

The  device,  circuit  and  process  designs  adopted  emphasize  future  growth  potential 
such  that  performance  characteristics  may  be  extended  to  much  higher  operating  frequen¬ 
cies  in  the  future. 


2.  THEORETICAL  ANALYSIS  OF  CCD  INPUT/OUTPUT  STRUCTURES 


2. 1  Purpose  and  Basic  Approach 

Considerable  effort  has  been  devoted  in  the  past  to  theoretical  analysis  of  the 
transport  of  charge  from  gate  to  gate  in  the  main  channel  of  a  CCD.  Relatively  little 
effort  has  been  applied  to  a  theoretical  evaluation  of  charge  flow  limitations  in  the  input 
and  output  structures  in  which  conversions  are  made  between  external  signals  and  charge 
in  the  CCD.  Some  of  the  results  of  the  charge  transfer  efficiency  studies  can  be  carried 
over,  but  there  are  significant  differences. 

Since  the  effects  of  charge  transfer  loss  accumulate  over  hundreds  or  even 
thousands  of  transfers,  residual  charges  in  the  order  of  10*  3,  10'^,  or  10*^  are  of  interest 
and  concern.  The  I/O  structures  typically  occur  only  once,  and  errors  do  not  accumulate. 
Charge  residues  affect  primarily  only  the  linearity  and  secondarily,  as  will  be  seen  later, 
the  noise  level  of  the  conversions.  Buried  channel  and  short  gate  structures  appropriate 
to  very  high  speed  CCDs  have  intrinsic  nonlinearities  of  a  few  percent  to  begin  with,  and 
thus  charge  residues  in  the  orders  10*1  and  10*2  are  of  interest  in  an  analysis  of  I/O 
performance. 

Another  significant  difference  is  that  the  transfer  of  charge  from  one  gate  to  the 
next  has  a  relatively  simple  formulation  with  boundary  conditions  that  permit  closed-form 
mathematical  expressions  to  be  found  for  special  cases.  The  input/output  structures  and 
processes  are  more  complex.  Numerical  analysis  is  therefore  the  primary  tool  for  solving 
the  equations,  although  the  results  of  the  numerical  analyses  may  suggest  simpler  and 
useful  ad  hoc  models. 

The  theoretical  problem  posed  by  the  real  CCD  I/O  structures  is  far  too  complex 
for  practical  solution  even  by  numerical  means.  A  high  speed  buried  channel  structure 
must  be  treated  in  at  least  two  dimensions,  and  the  time  required  to  design  and  run  such  a 
computer  program  would  be  very  great.  Even  then  the  relation  of  the  solutions  to  the 
laboratory  performance  of  devices  would  be  questionable  because  the  actual  dopant  distri¬ 
butions  in  multiple-implanted  and  oxidized  structures  is  not  well  known.  The  only  way  to 
get  reliable  specific  results  is  by  experimented  measurement.  The  purpose  of  any 
theoretical  analysis  should,  therefore,  be  to  gain  general  insight  into  the  physical  factors 
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that  limit  performance  and  on  their  general  dependence  on  the  dimensions  and  other 
design  variables  of  the  structure.  This  kind  of  information  is  useful  in  arriving  at 
intelligent  design  choices  and  in  suggesting  experimental  measurement  approaches. 

Since  the  real  problem  is  far  too  complex  to  analyze  practically,  we  will  develp  an 
admittedly  oversimplified  model  that  is  easy  to  formulate  as  a  computer  program  and 
requires  only  modest  computation  time.  The  results  will  not  be  quantitatively  accurate 
description  of  the  real  CCD,  but  they  will  show  trends  that  can  be  expected  to  be  reliable 
as  design  parameters  and  techniques  of  operation  are  varied. 

2.2  Formulation  of  Problem 

2.2.1  General  Formulation 

We  will  begin  with  a  more-or-less  complete  and  general  formulation  of  the 
charge  flow  problem  and  then  systematically  introduce  the  series  of  approximations  that 
leads  to  our  final  model. 

The  effects  of  electrostatic  fields  and  carrier  diffusion  can  be  handled 
together  using  the  Fermi  potential  function; 

°f  =  V  ’  Vthln(n/no)  (1) 

where  V  is  the  actual  electrostatic  potential,  Vt^  the  thermal  voltage  KT/qe,  n  the  den¬ 
sity  of  carriers  (electrons),  and  no  an  arbitrary  constant,  conveniently  chosen  in  some 
cases  to  be  the  intrinsic  carrier  concentration  nj  =  1.5  x  10^cm-3.  In  our  analysis  we  will 
ignore  minority  carriers,  since  in  all  cases  of  interest  n  is  much  greater  than  nj  and  p. 

From  the  Fermi  potential  we  derive  the  effective  electric  field  acting  on 

the  carriers. 


e  = 


-VO, 


(2) 


The  velocity  of  the  carriers  is  then  given  as; 


This  can  be  taken  as  the  definition  of  the  net  mobility  y.  In  practice,  the 
relation  between  v  and  e  is  not  linear,  especially  at  high  field  strengths  in  the  range  of 
1  V/ ym.  Figure  2-1  shows  experimental  values  of  the  net  mobility  ynet  and  incre¬ 
mental  mobility  yjncr  *or  high-quality  silicon. 


Figure  2-1  -  Electric  Field  versus  Mobility  in  Silicon 


v 
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(4a) 

(4b) 


Note  that  the  net  mobility  which  appears  in  our  equations,  falls  off  more 
gradually  than  the  incremental  mobility.  Moreover,  in  silicon  with  a  low-field  mobility  of 
500  cm^/v-sec,  as  is  more  common  in  CCDs,  the  curve  probably  would  look  like  the 
dashed  curve  in  Figure  2-1.  As  will  be  seen  in  Subsection  3.1  later,  when  we  have  very 
short  effective  gate  lengths  (  ~  1  ym),  we  must  begin  to  consider  the  effect  of  mobility 
reduction. 
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From  the  carrier  velocity  given  in  Equation  (3),  we  derive  the  flux  of  car¬ 
riers  F  as  follows: 

T  =  nv"  =  -np7  =  npyd.  (5) 

The  continuity  equation  expresses  the  fact  that  if  no  carriers  are  spontaneously  generated 
(a  reasonable  assumption  over  the  short  time  intervals  under  consideration  here),  the 
change  in  carrier  density  is  the  negative  divergence  of  the  flow;  that  is, 

-r—  =  -  V  .  F  =  -  V  .  (pnV9f) 

at 

=  -p  y.  ny<?.- 

Equation  (6b)  holds  when  p  is  constant,  as  we  will  assume  from  this  point  on.  Equation 
(6a)  is  fairly  generally  valid.  One  significant  assumption  not  yet  stated  explicitly  is  that 
electrodynamic  effects  are  ignored.  These  include  radiation  of  electromagnetic  fields 
from  accelerated  charges  and  magnetic  fields  induced  by  time-varying  electric  fields. 

2.2.2  Approximations  and  Simplifications 

2. 2. 2.1  Limitation  to  One  Dimension 

Our  first  major  simplification  is  to  treat  all  the  variables  as  a  func¬ 
tion  only  of  the  dimension  parallel  to  the  direction  of  charge  transport  in  the  CCD  chan¬ 
nel,  which  we  will  call  the  X-axis.  Distributions  of  carriers  into  the  depth  of  the  struc¬ 
ture  will  be  integrated  into  an  effective  area  concentration.  These  concentrations  will  be 
assumed  not  to  vary  across  the  width  of  the  device.  Since  the  width  will  typically  be 
50  um  while  the  gate  lengths  and  channel  depths  will  be  5  pm  or  less,  this  approxi¬ 
mation  is  quite  reasonable.  Carrier  motion  in  the  depth  direction  will  be  neglected. 


(6a) 

(6b) 


Equations  (2)  through  (6)  now  take  on  the  form: 


dn  _  JL 
3t  ~  ‘  Ox 


(7a) 


(7b) 


get: 


Inserting  the  definition  of  the  Fermi  potential  from  Equation  (1),  we 


3n  3 
3t  "  3x 


pVth 


(8) 


The  first  term  represents  field-driven  flow;  the  second  term  represents  thermal  diffusion 
(yVtf,  =  ykT/qe  =  Dtf,  by  the  Einstein  relationship). 

2.2.Z2  Local  Determination  of  Potential 

We  will  now  make  the  major  assumption  that  the  electrostatic  poten¬ 
tial  V(x,t)  at  a  given  point  x  depends  only  on  the  value  of  the  carrier  density  n(x,t)  at  that 
point.  This  is  certainly  not  always  true;  the  potential  in  general  depends  on  the  charge 
distribution  everywhere.  For  surface  channel  devices  with  the  signal  charge  separated 
from  the  gates  only  by  a  thin  oxide  of  perhaps  0.1  ym  and  with  gates  on  the  order  of  10 
ym  long,  the  approximation  is  very  good.  The  approximation  becomes  less  good  as  the 
channel  becomes  more  deeply  buried  (  1/2  to  1  ym  below  the  gates)  and  as  the  smallest 
geometrical  features  of  the  gates  shrink  towared  1  ym. 

The  mathematical  consequence  of  this  assumption  is  the  expression: 

av  .  9Vch  +  dV  on  (9) 

ax  ~  ax  an  ax 


where  Vch  is  the  channel  potential  when  n  =  0.  The  second  term  in  Equation  (9)  accounts 
for  the  change  in  potential  caused  by  the  local  carrier  density.  Stated  another  way,  the 
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first  term  represents  the  fringing  electric  field  from  the  gates  and  the  second  term  repre¬ 
sents  the  field  due  to  mutual  repulsion  of  the  charged  carriers.  Putting  Equation  (9)  into 
Equation  (8),  we  get: 


an 

aj 

I  3Vch 

av 

an  I 

at 

3x  | 

3x 

l 

-M- 

V  ,  -  n  - — 
th  3n 

► 

3x  j 

(10) 


The  term  u  -  n  j  represents  an  effective  diffusion  constant  in  which  the 

thermal  diffusion  constant  pVth  =  Dt^  is  augmented  by  |  iLX.  is  negative  j 


a  term 


proportional  to  the  local  carrier  concentration. 

2.2.2.3  Constant  Effective  Capacitance 


We  will  now  make  the  relatively  minor  assumption  of  a  constant 
effective  capacitance  in  the  relationship  between  V  and  n.  We  thus  write: 


ch 


q  n 
e 

C 


(11) 


and, 


3  V 
dn 


c 


(12) 


This  approximation  also  is  very  good  for  surface  channel  devices  where  the  total  effective 
capacitance  is  dominated  by  voltage-independent  oxide  capacitance.  In  a  buried  channel 
device,  the  effective  capacitance  of  a  nearly  empty  channel  is  significantly  smaller  than 
that  of  a  nearly  full  one.  One  would,  on  general  principles,  expect  improved  performance 
because  the  effective  diffusion  constant  at  lower  carrier  concentrations  would  be  larger, 
and  this  is  when  the  extra  diffusion  speed  helps  the  most. 

Substituting  Equation  (12)  into  Equation  (10)  we  get: 
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ax 

(13) 
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We  can  now  express  the  carrier  concentrations  in  terms  of  the  change  in  channel  potential 
by  making  the  substitution, 


V 

n 


or. 


n 


V 

n 


(14a) 


(14b) 


where  we  have  chosen  the  symbol  Vn  to  represent  the  carrier  density  expressed  as  an 
equivalent  voltage.  Equation  (13)  can  then  be  converted  to  a  particularly  aesthetic  form: 


9Vn 

JT  j  05) 

2.2.3  Computer  Formulation 
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We  must  now  convert  Equation  (15),  a  second  order  partial  differential 
equation,  to  a  form  suitable  for  numerical  analysis  on  a  digital  computer.  As  with  the 
preceding  treatment,  we  take  the  simplest  and  most  direct  approach.  We  define  the  val¬ 
ues  of  Vch  and  which  we  will  denote  as  VCH  and  VN  for  the  computer  variables,  on  an 
array  of  x  values  speced  by  DX.  The  quantity  inside  the  braces  of  Equation  (15)  is  the 
flux.  It  will  be  denoted  as  F  and  defined  on  an  array  of  x  values  midway  between  the 
points  on  which  VCH  and  VN  are  defined. 

The  derivatives  are  replaced  by  finite  differences,  and  the  factor  1/DX  is 
taken  outside  the  braces.  Thus  we  have: 


VTH  + 


-'*{(■ 

VNfl+I  )+VN  (I) 


(1+1/2)  .MU*  |  ymil  J  [VCH(I+D  -  VCH (I)] 

i 

[VN(I+1)-VN(1)]> 


(16) 


Note  that  we  have  used  for  VN  (1+1/2)  the  average  of  the  values  at  the  two  points  I  and 
1+1  at  which  VN  is  defined.  Also,  since  the  computer  functions  with  integral  indices  only, 
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we  actually  add  1/2  to  the  indices  of  F  in  Equation  (16)  and  use  F(I+1)  in  place  of  F(Ul/2). 
We  will  continue  to  use  the  latter  expression  here,  however,  because  it  makes  the  relation¬ 
ships  of  terms  clearer. 

The  solution  proceeds  as  a  two  step  iteration.  First,  the  current  values  of 
VN  and  VC  H  are  used  to  calculate  values  of  F.  Then  the  values  of  F  are  used  to  get  new 
values  for  VN  using  the  following  Fortran-type  equation: 


VN(I) 


v::ii )  - 


DT 

(DX)2 


F(It1/2)  -  F(I-  1/2) 


(17) 


where  DT  is  the  time  increment  for  each  step  of  the  iteration.  Note  that  the  equation 
has  a  scaling  factor  (DT/DX^).  This  means  that  as  long  as  the  original  differential  equa¬ 
tion  remains  valid,  reducing  all  distances  by  a  factor  of  two  results  in  an  identical  solu¬ 
tion  evolving  four  times  faster  in  time.  Relative  changes  in  the  mobility  also  can  be 
extracted  as  part  of  the  scale  factor  so  that  doubling  the  mobility  again  results  in  an 
identical  solution  but  one  that  evolves  two  time  faster  in  time. 

Boundary  conditions  must  be  provided  at  the  first  and  last  points.  The  spe¬ 
cific  choice  will  depend  on  the  problem,  but  there  are  two  common  types.  One  occurs 
when  the  problem  area  terminates  with  a  high,  steep  potential  barrier.  We  can  represent 
this  case  by  requiring  that  the  flux  F  at  the  boundary  be  zero.  The  second  common 
boundary  condition  occurs  when  the  problem  area  terminates  with  a  deep,  steep  potential 
well.  This  case  is  commonly  represented  by  requiring  that  VN  be  zero  at  the  boundary.  In 
our  formulation  of  the  problem,  the  boundaries  occur  between  values  of  VN  and  at  values 
of  F.  We  therefore  handle  this  case  by  setting  the  first  value  of  VN  outside  the  problem 
area  to  zero  and  writing  the  boundary  value  of  the  flux  following  Equation  (16)  as  shown 
in  Equation  (18). 

F(boundary)  *  MU*  [VTH-i-l/2  VN (boundary))]  *VN (boundary)  (18) 

2.3  Cases  Analyzed 

Actually  getting  the  computer  program  written,  debugged,  and  running  well  was 
more  difficult  than  expected.  In  particular,  although  the  iterative  procedure  is  quite 
simple  and  straightforward,  it  does  not  always  converge.  Considerable  time  was  required 
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to  learn  about  the  problems  that  could  develop  in  the  calculation  and  to  get  the  program 
to  trap  and/or  recover  from  the  errors.  With  DX  =  0.2  pm  and  MU  =  500  cm2/v-sec, 
time  steps  no  larger  than  1/2  psec  and  some  as  small  as  1  fsec  were  required  to  achieve 
convergence. 

The  first  case  that  we  chose  to  look  at  using  the  computer  program  was  the 
equilibration  phase  of  the  potential  equilibration  input  technique.  We  chose  this  case 
because  it  is  relatively  simple  (the  channel  potentials  and  boundary  conditions  are  con¬ 
stant  in  time)  and  because  some  published  results  indicate  that  this  input  method  begins 
to  degrade  in  the  50  to  100  MHz  range.  We  were  curious  to  see  if  the  computer  program 
would  predict  a  speed  improvement  using  a  1  pm  effective  gate  length  of  a  noval  transis¬ 
tor  structure  that  has  been  developed  at  Raytheon  called  an  ECMOJ  transistor.  The 
channel  length  of  this  device  is  defined  by  using  an  undercut  etch,  hence  the  acronym 
ECMOS  for  "etched  channel  MOS." 

The  structure  of  the  problem  is  shown  in  Figure  2-2.  The  source  of  charge  is  to  the 
right  of  the  problem  area,  and  the  first  stage  of  the  CCD  is  to  the  left.  The  problem  area 
itself  contains  two  regions  -  an  input  well  (IW)  and  a  reference  well  (RW).  The  channel 
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Figure  2-2  -  Potential  Equilibration  Potentials  Just  Prior  to  Equilibration 
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potential  in  RW  is  1  V,  and  the  channel  potential  in  IW  is  V1N  +  1  V,  where  V1N  is  a  varia¬ 
ble  input  voltage.  The  input  well  has  a  width  of  5  pm;  the  reference  region  has  a 
variable  width  XRW.  The  equilibration  phase  begins  with  the  entire  problem  area  filled 
with  charge  to  a  Fermi  potential  of  Ov.  The  left  boundary  condition  is  F=0;  the  right 
boundary  condition  is  a  flux  corresponding  to  zero  carrier  density  just  outside  the  problem 
area. 

2.3. 1  Test  Case 

With  VIN=0  we  have  the  classic  charge  transfer  efficiency  problem,  for 
which  some  closed-form  mathematical  solutions  are  available.  We  used  this  problem  as  a 
test  case  to  verify  that  the  numerical  iteration  procedure  was  giving  plausible  results. 
The  self -field-aided  diffusion  problem,  which  we  simulated  by  setting  VTH=0,  has  a  solu¬ 
tion  of  the  form  n  (x,t)  =  nx(x)nt(t),  where  the  carrier  distribution  has  a  shape  nx(x)  that 
does  not  change  in  time.  The  time  dependence  has  the  form: 

nt(tl  1  (19) 

which  can  be  rearranged  to  give: 


To  test  the  constancy  of  the  spatial  distribution,  we  first  laid  the  computer¬ 
generated  graphs  of  VN(I)/VN(l)  over  each  other.  They  were  indistinguishable  for  times 
of  0.25  nsec  or  more.  As  a  more  sensitive  check  we  compared  the  ratios  of  peak  charge 
value  VN(l)  to  the  total  integrated  charge.  The  values  for  t=0.5,  1  and  2  nsec  agreed  to 
four  decimal  places;  the  value  of  0.25  nsec  differed  by  0.1  percent.  The  results  shown  in 
Figure  2-3  confirm  the  theoretical  time  dependence  of  Equations  (19)  and  (20).  For 
0.25  nsec  <  t  <  2  nsec  the  maximum  deviation  in  VN(1)  was  0.0004  V.  The  time  to  make 
the  transition  from  the  square  distribution  at  t=0  to  the  time-invariant  distribution  is 
about  75  psec. 
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Figure  2-3  -  Theoretical  Time  Dependence  of  Self -Field- Aided  Diffusion 
2.3.2  Sample  Case 

The  time  evolution  of  the  carrier  distribution  for  a  sample  case  is  shown  in 
Figures  2-4  and  2-5.  Figure  2-4  shows  the  distribution  of  VN;  Figure  2-5  shows  the  distri¬ 
bution  of  the  Fermi  potential  with  the  constant  (see  Equation  (1))  chosen  to  give  (J>prO  at 
t=0.  The  length  of  the  reference  gate  is  5  ym  and  the  value  of  the  input  signal  2  V  in 
this  case.  One  sees  from  these  plots  that  the  distribution  in  the  input  well  is  always 
almost  constant.  Because  of  the  Vn  term  in  the  effective  diffusion  constant,  the  effec¬ 
tive  diffusion  in  the  input  well  is  much  larger  than  that  in  the  reference  region.  Conse¬ 
quently,  the  carriers  in  the  input  well  rearrange  rapidly  enough  to  maintain  an 
equipotential  distribution. 

As  a  result,  for  sufficiently  large  values  of  V1N  (a  volt  or  more,  for 
example),  the  IW  region  could  be  modeled  as  a  single  region  of  width  DX  with  a 
capacitance  that  is  larger  by  the  ratio  XIW/DX,  the  number  of  DX  intervals  in  the  input 
region. 
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Figure  2-4  -  Time  Evolution  of  Charge  Level  versus  Position 
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2-3.3  General  Results 


The  equilibration  phase  of  the  potential  equilibration  input  method  was 
studied  for  a  range  of  values  of  VIN  and  values  of  XRW  of  1,  2  and  5  ym.  Figure  2-6 
shows  the  results  for  a  1  pm  reference  gate  at  four  values  of  VIN,  namely  VIN=0,  1,  2, 
and  3  V,  for  times  up  to  2.5  nsec.  The  ordinate  shows  the  amount  of  charge,  expressed  in 
volts,  that  remains  in  the  input  well  in  excess  of  the  amount  that  would  be  there  after  an 
infinite  equilibration  time.  Not  included  in  the  figure  is  the  excess  charge  that  remains  in 
the  reference  gate  region.  When  the  first  stage  of  the  CCD  shift  register  turns  on,  all  of 
the  charge  in  IW  and  some  of  the  charge  in  RW  will  transfer  into  the  CCD. 


Figure  2-6  -  Residual  Charge  versus  Equilibration  Time 


The  excess  charge  in  the  input  region  above  the  fully  equilibrated  value  has 
two  effects  on  the  performance  of  the  input.  First,  it  changes  the  relationship  between 
the  input  signal  and  the  charge  injected  into  the  CCD.  It  always  has  the  effect  of  intro¬ 
ducing  what  is  commoly  called  a  "fat  zero."  This  means  that  even  when  the  input  signal  is 
set  to  a  value  to  give  zero  charge,  the  actual  charge  will  be  larger  than  zero.  This 
reduces  the  available  dynamic  range.  In  addition,  if  the  excess  charge  does  not  depend 
linearly  on  VIN,  a  nonlinearity  will  be  introduced  into  the  voltage-to-charge  conversion. 
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The  charge  versus  voltage  relationship  derived  from  Figure  2-6  is  shown  in 
Figure  2-7  for  equilibration  times  of  0.5  nsec,  1  nsec,  2.5  nsec,  and  infinity.  Only  the 
charge  in  the  input  well  is  included,  and  it  is  expressed  as  the  equivalent  input  voltage  for 
inifinite  time  equilibration.  A  "least  squares  straight  line"  fit  reveals  that  the  nonlinear¬ 
ity  introduced  even  at  0.5  nsec  is  only  1  or  2  percent,  not  a  really  s'gnificant  degradation 
when  one  considers  that  buried  channels  and  fringing  effects  int:oduce  nonlinearities  of 
this  magnitude  also. 

The  second  effect  of  the  excess  charge  on  performance  is  an  increase  in 

noise  level.  When  complete  equilibration  is  allowed  to  occur,  only  thermal  noise  is 

present.  Roughly  speaking,  this  results  in  fluctuations  on  the  order  of  (0.03  V)  in  the 
voltage  level  in  the  input  well.  For  input  signals  of  1  V  or  more,  the  excess  charge  has 
decayed  to  this  level  at  around  2-1/2  to  4  nsec. 

Figure  2-8  shows  results  similar  to  those  of  Figure  2-6  but  for  a  reference 
gate  length  of  2  pm.  Figure  2-9  shows  the  result  with  VIN=2  V  for  a  5  pm  reference 

gate.  The  curves  all  have  a  similar  appearance;  only  the  scale  of  the  time  axis  is 

changed.  Figure  2-10  shows  the  time  required  to  reach  a  specified  residual  charge  as  a 
function  of  reference  gate  length  for  an  input  signal  of  2  V.  Generally  speaking,  the 
speed  of  the  equilibration  process  increases  as  L'1. 

The  linear  dependence  on  L  may  be  surprising  at  first  in  view  of  the  general 
A  scaling  law  mentioned  in  Subsection  2.2.3.  This  can  be  clarified  as  follows.  If  we 
compare  two  structures,  one  with  a  1  pm  reference  gate  and  5  pm  input  well  and 
another  with  a  2  pm  reference  gate  and  10  pm  input  well,  the  latter  will  respond  four 
times  more  slowly.  However,  if  we  keep  the  input  well  at  5  pm  in  both  cases,  there  is 
relatively  only  half  as  much  charge  that  must  flow  across  the  reference  region.  This 
contributes  a  two-fold  speed  increase  so  that  the  net  dependence  on  length  is  L  rather 
than  iA 

We  hav-'  operated  our  current  high-speed  CCD,  which  has  a  7-1/2  um  refer¬ 
ence  gate  and  a  15  pm  input  well,  at  speeds  up  to  60  MHz  with  good  results.  The  scaled 
gat  shrinkage  from  7-1/2  to  5  pm  should  result  in  approximately  a  two-fold  speed 
increase  (this  change  scales  as  A.  In  addition,  provided  the  drift  region  that  is  part  of 
the  ECMOS  structure  does  no  interfere  with  the  operation  of  the  input,  the  effective 
shrinking  of  the  reference  gate  to  about  1  pm  should  provide  an  additional  factor  of  five 
speed  improvement.  The  limit  could,  therefore,  be  as  high  as  0.5 GHz. 
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Figure  2-7  -  Effective  Input  Voltage  versus  Input  Voltage 
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RESIDUAL  CHARGE  (V)  M  RESIOUAL  CHANGE  (VI 


Figure  2-10  -  Time  to  Reach  Residual  Charge  Levels  versus  Reference  Gate  Length 


2.3.4  Equilibration  with  the  ECMOS  Structure 

The  computer  model  was  extended  to  include  an  ECMOS  drift  region,  as 
shown  in  Figure  2-11.  The  purpose  of  the  calculation  was  to  gain  insight  into  the  effect 
of  the  drift  region  on  the  speed  of  equilibration.  The  dimension.,  used  as  shown  in  the 
figure. 

Qualitatively,  what  happens  is  that  first  the  charge  in  the  drift  region  flows 
back  into  the  source.  This  happens  relatively  rapidly  because  of  the  high  carrier  density. 
Once  the  potential  in  the  part  of  the  drift  region  that  abutts  the  reference  well  region  has 
fallen  below  the  potential  of  the  bottom  of  the  reference  well,  the  drift  region  ceases  to 
be  a  factor  in  the  equilibration  process.  Figure  2-12  shows  results  from  the  conventional 
short  channel  reference  gate  and  for  an  ECMOS  gate  with  drift  region  potential  well 
depths  of  lv  and  3v.  There  is  an  initial  time  delay  with  the  ECMOS  gate  while  the  charge 
in  the  drift  region  subsides.  Apparently,  the  charge  distribution  in  the  input  and  refer¬ 
ence  wells  is  such,  however,  that  the  equilibration  process  quickly  catches  up.  For  Vp>  = 
lv,  the  results  are  essentially  identical  after  about  1.5  ns;  for  VD  =  3v,  less  than  half  a 
nanosecond  is  needed  to  catch  up.  Clearly,  an  ECMOS  input  gate  is  just  as  fast  and  effec¬ 
tive  as  a  true  short  gate. 


Figure  2-11  -  Structures  used  for  Modes  cf  Conventional  Type  Input  (A) 
and  ECMOS  Type  Input  (B) 
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Figure  2-12  -  Equilibration  Dynamics  for  ECMOS  input  Structures  Compared  to 
that  for  Conventional  Short-Channel  Structures 


2.3.5  Election  Velocity  Saturation 

The  calculations  presented  thus  far  have  assumed  a  constant  election 
mobility  regardless  of  electric  field  strength.  The  model  was  modified  to  check  the 
effect  of  election  velocity  saturation  by  replacing  the  carrier  velocity  by  a  constant  value 
VSAT  whenever  yE  was  larger  than  V5i^j.  In  general,  no  effect  was  found  when  the 
actual  experimental  value  of  saturation  velocity,  1  O'7  cm/s,  was  used.  To  verify  that  the 
model  did  include  the  effect,  some  runs  were  made  using  a  value  of  velocity  only  one  fifth 
of  the  actual  velocity.  Figure  2-13  shows  a  typical  result.  Early  in  the  equilibration,  a 
large  charge  gradient  -  and  hence  electric  field  -  exists  in  the  reference  well.  This 
gradient  relaxes  away  quickly,  and  soon  the  carrier  velocity  is  no  longer  saturated.  With 
a  velocity  saturation  value  of  10 7  cm/s,  a  voltage  difference  of  more  than  Iv  must  occur 
across  the  1  ym  reference  well  before  saturation  effects  will  be  significant.  For  a  one 
volt  overdrive  during  the  fill  phase,  this  condition  will  never  occur. 
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Figure  2-13  -  A  Significant  Velocity  Saturation  Effect 


2.4  Future  Theoretical  Work 

The  theoretical  work  carried  out  in  the  contract  has  consisted  of  the  development 
of  the  model  and  the  structure  of  the  computer  program,  and  the  computation  of  the 
performance  of  the  potential  equilibration  input  technique  during  the  equilibration  phase. 
Future  work  could  examine  additional  input/output  processes  and  the  effects  of  the 
approximations  used  in  the  model. 

2.4.1  Improved  ECMOS  Model 

The  model  of  the  ECMOS  input  in  Figure  2-11  was  not  complete.  It  included 
the  drift  region,  which  could  have  been  expected  possibly  to  have  an  effect  on  perform¬ 
ance,  but  it  neglected  the  fact  that  the  input  well  would  have  a  part  which  is  a  floating 
diffusion.  Figure  2-14  shows  the  structure  which  can  be  compared  to  Figure  2-1  i.  One 
would  not  expect  the  floating  diffusion  region  to  have  any  significant  effect  on 
performance. 


2-21 


12816 


CCD 

SHIFT. 

REGISTER 


FLOATING 

OlFFUSION 


INPUT  WELL 


ECMOS 

DRIFT  REGION 


H 


ECMOS  ACTIVE  REGION 


INPUT 

DIFFUSION 


Figure  2-14  -  ECMOS  Equilibration  Structure  Potential  versus  Position 
2-4.2  Injection  Model 

The  analysis  described  above  began  with  an  examination  of  the  equilibration 
phase  of  the  potential  equilibration  input  because  this  was  expected  to  be  the  performance- 
limiting  phase.  For  completeness,  the  injection  phase  should  also  be  investigated.  The 
same  computer  program  can  be  used  with  only  the  initial  conditions  (no  charge  present) 
and  the  boundary  condition  at  the  input  diode  side  changed.  Both  the  standard  MOS  and 
the  EC  MOS  structure  will  be  considered. 

2,4.3  Diode  Cutoff  Input  Model 

The  diode  cutoff  technique  is  expected  to  be  much  faster  than  the  potential 
equilibration  technique  because  the  charge  is  always  field-driven.  The  injection  phase  is 
similar  to  the  potential  equilibration  injection  phase.  The  cutoff  phase,  however,  involves 
a  time-dependent  change  is  the  reference  gate  level  from  at  or  below  the  input  well  level 
to  above  the  signal  level.  The  exact  way  in  which  the  cutoff  occurs  can  be  expected  to 
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influence  the  results.  Figure  2-15  shows  schematically  the  kinds  of  initial  intermediate, 
and  final  conditions  expected  with  the  standard  MOD  structure.  The  ECMOS  situation  is 
similar  but  more  complex. 


Figure  2-15  -  Diode  Input  Cutoff  Method  Operating  Conditions 


2.4.4  Output  Structures 

The  floating-diffusion  output  structure  can  be  operated  in  two  ways  which 
are  similar  to  the  input  techniques.  The  ECMOS  output  structure  is  essentially  identical 
to  the  input  structure;  the  conventional  MOS  output  structure  has,  in  addition,  the  float¬ 
ing  diffusion,  which  is  not  present  in  the  conventional  MOS  input  structure.  Conse¬ 
quently,  the  analysis  of  the  output  performance  should  be  straightforward. 

2.4.5  Effects  of  Approximations  in  the  Formulation 

A  number  of  rather  extreme  approximations  were  used  in  Subsection  2.2.2 
to  derive  the  formulation  that  will  be  used  for  the  bulk  of  the  theoretical  analysis.  Some 
insight  into  the  effects  of  these  approximations  can  be  achieved  by  treating  some  cases 
more  carefully. 
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The  simplest  improvement  to  make  is  to  introduce  gradations  in  the  values 
of  the  zero-charge  channel  potentials  to  represent  more  accurately  the  effects  of  fringing 
fields  from  the  gates  to  the  buried  channel. 

The  effect  of  carrier  velocity  saturation  and  voltage-dependent  capacitance 
can  also  be  examined  in  a  straightforward  manner.  The  nonlocal  dependence  of  electro¬ 
static  potential  on  charge  density  can  be  simulated  by  developing  the  next  order  of 
approximation  in  which  the  potential  depends  not  only  on  the  local  charge  density  but  also 
on  the  charge  density  in  the  neighboring  intervals. 

All  of  these  changes  will  add  considerable  time  to  the  core  computation  in 
the  innermost  loop.  Computer  time  may  become  too  great  to  permit  exhaustive  analyses 
to  be  carried  out  with  the  more  exact  formulation,  but  any  significant  trends  should  be 
revealed  from  a  few  key  cases. 


3.  CCD/ECMOS  COMPATIBLE  TECHNOLOGY 


During  this  phase  of  the  contract,  efforts  were  concentrateo  on  working  out  design 
and  processing  details  for  the  high  speed  CCD  as  well  as  the  ECMOS  devices  to  be  used 
for  the  CCD  input,  output,  output  amplifier,  and  clock  driving  devices. 

The  design  and  fabrication  of  the  high  speed  CCD  delay  line  closely  followed  the 
methods  successfully  used  for  the  Raytheon  60  MHz  CCD.  The  frequency  limitations  of 
this  device  at  60  MHz  are  primarily  due  to  the  output  amplifier  response.  The  frequency 
response  limit  imposed  by  the  charge  transfer  properties  of  the  CCD  is  substantially 
higher. 

The  above  conclusion  was  derived  from  the  fact  that  the  output  pulse  responses  for 
the  64  stage  and  128  stage  devices  looked  identical  when  clocked  at  60  MHz.  The  excel¬ 
lent  high  frequency  performance  of  these  CCDs  was  obtained  in  spite  of  the  fact  that  the 
gate  length  was  7.5  pm  for  the  transfer  and  storage  ceils,  and  the  fact  that  the  devices 
were  of  the  two-phase  type.  In  addition,  all  gate  electrodes  were  made  of  polysilicon, 
with  electrode  overlaps  of  1.25  pm.  For  the  purpose  of  this  contract,  the  channel  width 
of  the  CCD  was  reduced  from  225  pm  to  a  value  in  the  range  50-100  pm,  and  cell 
electrode  lengths  lowered  to  5  pm.  The  gate  electrode  overlap  of  1.25  pm  was 
retained.  These  changes  were  implemented  to  improve  high  frequency  transfer  efficiency 
and  reduce  the  clock  capacitance  to  be  driven.  A  further  modification  was  the  use  of 
poly  Si- A  l  or  TiW-AI  instead  of  polysilicon  for  the  second  layer  of  gate  electrodes 
(storage  cell  electrodes).  The  reduced  cell  electrode  lengths  enable  the  straight  256  stage 
device  to  be  built  with  33  percent  shorter  overall  length. 

In  the  past,  the  60  MHz  CCDs  have  been  made  using  an  ion  implanted  guard  band 
to  confine  the  charge  to  the  CCD  area  in  the  direction  perpendicular  to  its  propagation 
and  to  isolate  the  transistors  of  the  amplifier  from  each  other.  This  isolation  method  uses 
up  more  space,  requires  a  separate  registration  step,  and  leads  to  higher  oxide  steps  than 
the  so  called  coplanar  or  semirox  isolation  methods. 

In  the  coplanar  method,  the  device  to  be  isolated  is  first  protected  by  an  ozidation 
mask  made  up  of  Si02  and  SijNtf.  The  unmasked  silicon  substrate  is  etched  to  a  depth 
about  one  half  the  desired  field  oxide  thickness.  Boron  is  then  implanted  into  the  surface 
of  the  high  resistivity  P-type  silicon  to  prevent  inversion,  and  field  oxide  is  then  grown  by 
thermal  oxidation.  After  this,  the  masking  layer  is  stripped  and  the  device  is  formed. 
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The  semirox  isolation  method  is  a  modification  of  the  coplanar  technique  which 
omits  the  etching  of  the  substrate.  As  a  result,  the  fieid  oxide  forms  a  step  which  is 
approximately  one  half  of  its  oxide  thickness  due  to  the  silicon  consumed  during  the 
thermal  oxidation. 

The  coplanar  and  semirox  isolation  methods  offer  the  important  advantage  that 
they  permit  direct  contact  of  field  regions  to  heavily  doped  regions  such  as  sources  and 
drains  of  field  effect  transistors  without  causing  low  reverse  breakdown  voltages.  This  is 
not  possible  with  conventional  guard  band  isolation  since  guard  bands  are  heavily  doped 
and  adequate  spacing  must  be  provided  between  them  and  the  source  and  drain  regions. 

While  ECMOS  transistors  are  readily  compatible  with  the  coplanar  or  semirox  isola¬ 
tion,  buried  channel  CCDs  using  Raytheons  existing  process  or  not.  While  surface  channel 
CCDs  have  been  fabricated  at  Raytheon  with  the  coplanar  isolation  which  perform  as  well 
as  guard  banded  devices,  we  have  not  achieved  similar  results  with  buried  channel  CCDs. 
It  is  believed  the  phosphorus  implant  used  to  create  the  buried  channel  penetrates  the 
fringe  of  the  isolation  region  and  overcompensates  its  boron  doping.  As  a  consequence, 
spurious  channels  have  been  formed  at  the  boundaries  of  the  CCD  buried  channels  short¬ 
ing  them  out. 

Although  this  situation  could  be  corrected  by  increasing  the  boron  dose  in  the  field, 
this  would  also  reduce  the  breakdown  voltage  of  the  source  and  drain  regions.  For  this 
reason,  the  new  high  speed  CCD  test  pattern  will  utilize  coplanar  isolation  for  I/O  and 
transistor  structures  and  guard  band  isolation  for  the  buried  channel  CCDs. 

3.1  ECMOS  Transistor  Device  Development 

Although  the  work  described  below  was  not  supported  by  contract  funds,  this  sec¬ 
tion  is  included  since  ECMOS  transistors  will  be  used  in  the  high  speed  CCD/ECMOS 
technology  under  development.  The  ECMOS  transistor  offers  the  main  features  of  the 
well  known  D-MOS  and  V-MOS  transistors,  namely  a  short  active  channel  length  (of  the 
order  of  1  pm)  which  does  not  require  precision  lithography  for  its  fabrication,  and  a 
relatively  high  drain  to  source  breakdown  voltage.  The  main  advantages  of  ECMOS  rela¬ 
tive  to  this  contract  reside  in  its  excellent  threshold  control,  small  gate  capacitance,  and 
the  compatibility  of  the  fabrication  procedure  with  that  of  CCDs. 

The  gate  capacitance  of  ECMOS  transistors  can  be  relatively  small  since  the  por¬ 
tion  of  the  gate  covering  the  drift  region  can  be  separated  from  the  substrate  by  an  oxide 
which  is  substantially  thicker  than  the  gate  oxide. 
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As  with  D-MOS  and  V-MOS  transistors,  a  high  dram  to  source  breakdown  voltage  is 
achieved  by  employing  a  drift  region  in  series  with  the  short  channel  control  regions.  The 
purpose  of  the  drift  region  is  to  accommodate  the  voltage  difference  between  the  drain 
voltage  and  the  voltage  which  can  be  supported  across  the  active  short  channel  region 
without  punching  through.  The  punch-through  voltage  increases  with  the  square  root  of 
the  p-doping  concentration  in  the  short  channel  region  and  thus  increases  with  the  thresh¬ 
old  voltage  for  a  given  gate  oxide  thickness. 

In  order  to  obtain  optimum  performance  from  such  a  transistor  the  electric  field  in 
the  short  channel  region  should  approach  10^  V/cm.  At  this  electric  field  strength  the 
carrier  velocity  starts  to  saturate  and  close  to  the  minimum  possible  transit  time  is 
achieved.  A  higher  field  strength  is  not  advisable  since  it  will  not  greatly  increase  carrier 
velocity  and  will  aggravate  short  channel  effects  (drain-voltage  dependence  of  the  thresh¬ 
old  voltage,  and/or  punch-through). 

The  family  of  drain  characteristics  with  gate  voltage  as  a  parameter  can  be  used 
to  determine  whether  or  not  a  given  device  attains  optimum  carrier  velocity.  As  the  drift 
velocity  saturates  the  transconductance  reaches  its  upper  limit  and  becomes  independent 
of  the  gate  voltage.  This  can  be  shown  as  follows  in  Equations  (21)  and  (22). 
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where 

W  =  channel  width 

L  =  channel  length 

vs  =  saturated  drift  velocity 

T  =  transit  time 

Q  =  charge  induced  by  Vq  -Vt 

Vq  =  gate  voltage 

Vt  =  threshold  voltage 

t  =  thickness  of  gate  oxide 
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Note  in  Equation  (22)  the  interesting  fact  that  at  velocity  saturation  the  drain 
current  and  the  transconductance  ao  not  depend  on  the  channel  length,  L.  However,  a 
channel  length  shorter  than  necessary  to  achieve  velocity  saturation  can  still  be 
important  for  high  frequency  operation  since  it  affects  the  transmit  time 

^  ,  the  gate  capacitance  ^ox  =  ~ t  '°  (W) 

and  the  on-resistance  of  such  a  device. 

The  parameters  for  an  ECMOS  test  structure  can  be  used  to  determine  whether 
they  operate  at  close  to  the  carrier  velocity  limit.  Using  t  =  7  x  10'^cm,  W  =  7.5  x  10*^ 
cm  and  <gm^max  -  1-8  mA/V  in  Equation  (22)  we  obtain  a  vs  of  4.75  x  10*>  cm/sec.  This  is 
approximately  80  percent  of  the  saturated  velocity  (6  x  10^  cm/sec)  available  in  bulk 
silicon. 

The  following  requirements  apply  concerning  the  drift  region: 

1)  It  must  conduct  sufficiently  to  permit  attainment  of  an  electric  field 
Es  =  1  x  10^  V/cm  in  the  active  channel  region 

2)  It  must  be  capable  of  carrying  the  current  delivered  by  the  short 
channel  region  under  optimum  conditions 


Assuming  an  average  drift  velocity,  v^ 
implies  Equation  (23a)  which  can  be  as  (23b). 


for  the  drift  region,  the  second  condition 


(23a) 


(23b) 


where 


t^  =  oxide  thickness  over  the  drift  region 
Vtcj  =  threshold  voltage  in  the  drift  region 
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For  practical  case  one  may  assume 


v  ,  «  0.  1  v  i  and  —  ss  3 

d  3  ’  : 

therefore 

VG‘Vtd  5830  (VG‘Vt)  {2U) 

Consequently,  if  a  gate  voltage  exceeding  Vt  by  10  V  is  desired  (corresponding  to  ID/W  = 
3A/'cm)  the  threshold  voltage  for  the  drift  region  must  be  negative  and  Vtd  *  -300  V. 
This  corresponds  to  the  charge  density  shown  in  Equation  (25)  which  corresponds  to  a 
density  of  electronics  charge  of  N  =  2 .4  x  10^/cm2. 

Q/cm2  =  ^-2*  V  ,  =  *--3'  85  *  l0, -  (300)  =  3.  3  x  10*6ccul/cm2  (25) 

d  za  2,  8  x  10* 3 

The  charge  density  calculated  in  Equation  (25)  can  be  located  in  the  oxide  as  a 
positive  fixed  charge,  or,  more  conveniently,  can  have  the  form  of  a  donor  implant  into 
the  surface  of  the  silicon.  For  both  the  value  of  N  is  rather  high  to  achieve  in  practice. 
More  easily  attainable  conditions  require  that  an  N  of  approximately  0.6  x  1012/cm2  be 
used  which  will  raise  vd  to  approximately  0.4  vs. 

Assuming  a  phosphorus  implant  is  used  to  obtain  the  above  conditions  a  dose  of  6  x 
10^2/cm2  and  a  junction  depth  of  7  x  10"^  cm,  the  dopant  concentration  will  be  in  the 
order  of  8.6  x  10  ^/cm  \ 

A  boron  concentration  in  excess  of  this  value  must  then  be  added  to  provide 
threshold  control  in  the  short -channel  region,  and  to  electrically  isolate  the  source  of  the 
ECMOS  transistor  as  shown  in  Figure  3-3. 

Figures  3-1  through  3-3  show  three  different  ECMOS  transistor  structures  under 
investigation.  Structure  (Figure  3-1)  will  be  referred  to  as  type  A.  It  relies  on  the  fixed 
positive  charge  contained  in  the  oxide  covering  the  drift  region  Lq  to  invert  the  drift 
region  more  strongly  than  the  channel  region  at  low  gate  voltages.  The  structure  shown 
in  Figure  3-2  as  type  B  has  a  phosphorus  implant  in  most  of  the  drift  region. 

Structure  (Figure  3-3,  type  C)  is  the  most  suitable  one  of  the  three  for  output 
amplifier  applications.  In  this  case  the  phosphorus  doping  is  implanted  during  the  CCD 
fabrication  sequence  and  before  the  device  areas  are  formed.  It  does  not  require  accu¬ 
rate  masking  since  the  excess  implant  is  automatically  removed  when  the  silicon  surface 
is  etched  prior  to  the  formation  of  the  coplanar  field  oxide. 
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Figure  3-1  -  ECMOS  Transistor  Suitable  for  Small  Gate  Voltage  Applications 


Note  that  the  boron  threshold  implant  serves  two  functions  in  the  type  C  structure. 
It  serves  to  adjust  the  threshold  for  the  short  channel  region,  and  also  to  diode-isolate  the 
source  from  the  n-type  layer.  To  serve  both  functions  the  boron  must  therefore  be  driven 
to  a  depth  exceeding  the  shallow  arsenic-doped  source  region.  Typical  S/D  families  and 
surface  potential  vs  gate  voltage  characteristics  are  shown  in  Figures  3-4  through  3-6  for 
the  ECMOS  structures  shown  in  Figures  3-1  and  3-3.  In  all  cases,  the  length  of  the  short 
channel  section  was  L  =  1  ym,  while  the  width  was  W  =  7.5  ym.  For  aach  type  of 
device,  the  effect  of  having  a  drift  region  of  either  Lq  =  4  ym  or  Lq  =  6.5  ym  is  also 
indicated. 

As  may  be  expected  from  the  previous  discussion,  Figures  3-4  and  3-5  show  a 
dramatic  change  in  device  characteristics  with  the  length  of  the  drift  region.  When  the 
phosphorus  dose  is  relatively  small  (1  x  10 12/cm 2)  an(j  Lq  =  6.5  ym  the  electric  field  at 
the  short  channel  region  is  too  small  to  support  a  saturated  carrier  velocity.  In  contrast 
for  a  phosphorus  dose  of  3  x  lO^/cm^,  there  is  much  less  dependence  on  Lq  since  the 
drift  region  impedance  is  much  less  than  that  of  the  short  channel  region. 

3.2  CCD/ECMOS  Wafer  Fabrication  Process  Development 

The  major  features  of  the  CCD  fabrication  method  to  be  utilized  can  be  summa¬ 
rized  as  follows: 

1)  Ion  implantation  will  be  used  for  most  doping  steps  including  isola¬ 
tion,  gettering,  buried  channel,  threshold  control,  and  source/drain 
because  of  the  wide  degree  of  accuracy  and  reproducibility  it 
affords,  as  well  as  the  cleanliness  it  achieves  by  doping  though  an 
oxide  and  the  convenience  of  using  photo-resist  as  an  implantation 
mask. 

2)  The  potential  profile  for  2  phase  buried  channel  CCDs  will  be  created 
by  the  use  of  two  phosphorus  implants  avoiding  the  generally  used 
and  difficult  to  control  compensation  method  using  phosphorus  and 
boron  dopants.  In  addition,  the  amount  of  dopant  in  the  substrate 
will  be  minimized. 

3)  Self  aligned  gates  will  be  created  by  doping  source  and  drain  regions 
while  they  are  covered  by  the  gate  oxide.  Thus,  undercutting  of  the 
gate  oxide  under  the  gate  electrodes  will  be  avoided. 
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Figure  3-5  -  Source/Drain  Families  of  ECMOS  Devices  (Type  C) 
with  Lj>  =  4  ym  and  6.5  ym 
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Figure  3-6  -  Gate  Voltage  versus  Surface  Potential  with  Nt2  as  a  Parameter 


4)  The  polysilicon  layer  will  be  doped  by  phosphorus  diffusion  since  this 
results  in  a  very  low  sheet  resistance  (20  £22  for  5000  A  of 
thickness). 

5)  ECMOS  technology  will  be  used  for  CCD  input  and  output  stages  to 
achieve  the  performance  of  micron  and  submicron  active  channels 
without  loss  of  charge  handling  capacity. 

6)  The  substrate  will  be  gettered  by  a  BF2  implant  on  the  back,  which 
also  will  provide  a  good  back  (substrate)  contact. 

7)  Polysilicon-aluminum  or  TiW-Al  metallization  will  be  used  to  provide 
uniform  step  coverage  and  reliable  contact  to  very  shallow  junctions. 

8)  Mixed  coplanar  and  guard  band  isolation  will  be  used  to  optimize  the 
performance  and  ease  of  fabrication  of  the  CCD  delay  line  and  I/O 
regions  ind  peripheral  transistor  circuitry. 

9)  Ultraclean  wafer  fabrication  processing  will  be  used,  including  the 
use  of  polysilicon  furnace  tubes  cleaned  with  HC1  and  tight  control 
over  the  quality  of  doping  and  gas  supplies. 
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Figure  3-9  -  Steps  5  and  6  of  CCD/ECMOS  Compatible  Process 
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Figure  3-10  -  Steps  8  through  10  of  CCD/ECMOS  Compatible  Process 
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Figure  3-11  -  Steps  11  through  14  of  CCD/ECMOS  Compatible  Process 
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Figure  3-12  -  Steps  15  through  IS  of  CCD/ECMOS  Compatible  Process 


Figure  *-13  -  Steps  19  through  25  of  CCD/ECMOS  Compatible  Process 


4.  HIGH  SPEED  CCD/ECMOS  MASK  DESIGN 

The  proposed  advanced  high  speed  CCD  technology  includes  simultaneously  fabri¬ 
cated  high  speed  CCD  delay  lines,  higher  speed  input/output  structures  and  high  fre¬ 
quency  MOS  transistors  for  both  low  voltage  amplifier  and  high  voltage  clock  driving 
applications. 

A  new  mask  series  was  designed  to  evaluate  various  approaches  to  such  a  technol¬ 
ogy.  The  mask  series  incorporates  input/output  test  structures  for  the  evaluation  of  con¬ 
ventional  high  performance  MOS  and  EC  MOS  gate  structures  with  a  variety  of  channel 
lengths  and  widths. 

Figure  4-1  shows  the  functional  representation  of  an  EC  MOS  input/output  test 
structure.  Note  that  it  amounts  to  a  high  speed  CCD  with  only  a  single  storage  well. 

This  pattern  will  permit  the  characterization  of  EC  MOS  I/O  structures  for  a  high  speed 

CCD  without  the  signal  degradation  caused  by  multiple  transfer  and  storage  stages.  f 


Figure  4-2  is  a  photograph  of  the  above  function  using  channel  widths  of  12.5,  25, 
50,  100  and  200  Urn.  The  input  and  output  electrodes  used  are  of  the  ECMOS  type  and 
all  inputs  and  outputs  are  brought  to  external  pads  so  that  they  can  be  evaluated  under  a 
variety  of  bias  and  signal  conditions.  The  purpose  of  this  group  of  test  patterns  is  to  pro¬ 
vide  comparison  data  for  five  channel  widths  with  respect  to  harmonic  distortion  noise, 
and  transfer  efficiency  at  both  low  and  high  operating  speeds. 

Figure  4-3  shows  the  functional  representation  of  a  conventional  polysilicon  gate 
input/output  test  structure. 

Figure  4-4  is  a  photograph  of  the  above  function  using  electrode  lengths  of  5,  10, 
and  40  ym  and  a  200  ym  channel  width. 

In  addition  to  the  test  cells  shown  in  Figures  4-2  and  4-4,  the  chip  will  contain  one 
256  stage  CCD  using  ECMOS  devices  for  I/O  electrodes  and  output  amplifier  devices  and 
a  second  256  stage  CCD  using  identical  5  ym  transfer  and  storage  electrodes  but  with 
more  conventional  I/O  and  amplifier  structures.  Performance  of  a  given  CCD  with  con¬ 
ventional  polysilicon  input /output  structures  will  therefore  be  directly  comparable  with 
the  performance  of  the  same  CCD  using  ECMOS  I/O  and  output  amplifier  structures. 
Both  256  stage  CCDs  will  use  50  ym  wide  channels  and  have  5  ym  polysilicon  transfer 
wells  and  5  ym  wide  aluminum  storage  wells. 

It  is  also  desirable  to  evaluate  CCD  performance  using  storage  and  transfer  wells 
that  are  not  equal  in  length.  In  order  to  evaluate  this  option,  three  additional  64  stage 
CCD  delay  lines  with  varying  lengths  for  transfer  and  storage  wells  will  be  provided  for 
this  purpose  in  the  mask  series. 

The  chip  was  divided  into  the  following  four  sections: 

a)  Section  A  contains  a  256  stage  CCD  with  ECMOS  drivers,  S/H,  I/O 
and  amplifier  transistor  structures. 

b)  Section  B  contains  an  identical  256  stage  CCD,  but  with  drivers,  S/H, 
I/O  and  amplifier  circuitry  constructed  from  conventional  self- 
aligned  polygate  NMOS  transistors. 

c)  Section  C  contains  input  and  output  CCD  test  structures  of  varying 
widths  and  lengths  containing  both  ECMOS  and  conventional  transis¬ 
tors.  Section  C  also  contains  ring  oscillators  and  individual 
transistors. 
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Figure  4-2  -  ECMOS  Test  CCD's  W=100,  200  u  m 


Figure  4-3  -  Functional  Representation  of  a  Conventional  Polysilicon 
Gate  Input/Output  Test  Structure 
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Figure  4-4  -  Conventional  NMOS 


d)  Section  D  contains  a  128  stage  CCD  without  drivers  and  with  conven¬ 
tional  devices  for  input/output  structures,  S/H  and  output  amplifier. 
Section  D  also  contains  ring  oscillators  and  individual  components. 

Each  of  these  four  sections  are  bounded  with  a  4  mil  grid  and  the  sections  can  be 
scribed  and  separated  individually.  The  overall  chip  size  is  177.6  x  241.8  mils. 

The  combined  ECMOS,  CCD  and  conventional  self -aligned  polyprocess  requires  12 
mask  layers  through  aluminum  metallization.  Additional  options  for  evaluation  of  an  N+ 
implant  in  the  ECMOS  drift  regions  of  the  input  and  output  transistor  structures  are  pro¬ 
vided.  Metal  mask  options  for  the  evaluation  of  drivers  and  amplifiers  independent  of  the 
CCD  are  also  provided.  Finally,  an  option  for  mask  defined  short  channel  (1  1/4  p) 
polygate  structures  is  provided.  This  requires  a  total  of  16  reticles  to  be  generated  and 
stepped  into  the  mask  array.  Evaluation  of  the  short  channel  (1  1/4  p)  conventional 
devices  will  require  that  wafers  are  processed  separately  for  that  purpose.  That  is,  all 
devices  on  the  self-aligned  polysixth  mask  will  be  short  channel.  The  N  implant  under 
the  I/O  structures  and  the  individual  device  metal  options  will  be  provided  by  stepping 
alternate  rows  of  options  into  the  final  masks. 

There  are  a  total  of  382  bonding  pads  on  the  chip,  arranged  in  10  rows  of  4  mil  pads 
on  6  mil  centers.  Each  of  the  256  stage  CCD  arrays  are  arranged  such  that  the  bonding 
pads  will  accommodate  useage  of  the  same  probe  card.  The  remaining  6  rows  of  bonding 
pads  comprise  a  random  arrangement  of  input/output  and  test  structures  that  will  require 
additional  probe  cards.  It  is  expected  that  a  single  curve  tracer  or  DC  probing  evaluation 
would  be  more  than  adequate  to  determine  whether  or  not  the  test  devices  should  be 
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assembled.  Evaluation  of  high  frequency  performance  will  require  the  assembly  of  test 
devices  into  a  package  arrangement  capable  of  handling  the  power  and  frequency 
requirements. 

In  addition  to  the  CCD  arrays,  ring  oscillators,  test  transistors,  resistors  and 
capacitors  were  added  to  the  chip.  The  ring  oscillators  consist  of  two  ECMOS  and  two 
standard  polygate  structures  of  19  stages  each.  An  additional  two  stages  are  added  to 
each  oscillator  to  buffer  the  oscillator  stages  from  the  output  terminal.  Both  of  the 
ECMOS  ring  oscillators  and  the  standard  ring  oscillators  are  made  with  depletion  loads 
and  with  a  separate  terminal  that  can  be  used  to  connect  all  the  gates  of  the  load  devices 
to  a  Vjjb  bias  terminal.  These  ring  oscillators  can  be  used  to  provide  an  indication  of  the 
gate  delays  between  ECMOS  and  standard  polysilicon  gates. 

The  capacitors  are  constructed  from  the  various  implants  to  provide  information 
on  the  oxides  of  the  process  steps.  Four  point  resistors  are  also  provided  to  monitor 
the  sheet  resistance  of  the  implants.  Finally,  test  transistor  structures  are  provided 
for  evaluating  the  performance  of  the  ECMOS  and  polygate  devices. 
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5.  RESULTS 


5. 1  Introduction 

Unfortunately,  two  unexplained  processing  problems  prevented  the  demonstration 
operable  devices. 

The  first  problem  was  the  disappearance  of  the  nitride  layer  used  for  the  con¬ 
trolled  undercut  etching  that  forms  the  ECMOS  transistors.  At  no  time  were  the  wafers 
knowingly  placed  in  a  nitride  etch,  but  at  the  critical  processing  step  it  was  found  that 
the  nitride  layer  had  vanished.  The  conventional  type  MOS  structures,  however,  were 
expected  to  be  usable. 

The  second  problem  was  an  oxide  problem  never  seen  before.  Transistors  whose 
threshold  voltages  were  normally  within  a  volt  of  zero  exhibited  threshold  voltages  of 
about  -5  V  or  -10  V.  Batches  of  wafers  processed  together  showed  some  wafers  with  one 
value  and  some  with  the  other.  It  would  appear  that  the  starting  material  was  defective 
except  that  the  buried  channel  regions  had  the  threshold  characteristics  expected.  These 
wafers  also  exhibited  very  high  dark  currents  and  other  peculiarities. 

It  appeared  that  the  experimental  part  of  this  program  was  going  to  yield  no  useful 
information.  The  one  chance  to  learn  something  useful  lay  in  the  set  of  three  single-stage 
test  CCD's  with  different  size  input  gates  (see  Table  5-1).  These  devices  were  included 

TABLE  5-1 

CCD’s  WITH  DIFFERENT  SIZE  INPUT  GATES 


Device  Type 

wmmmm 

ZGl 

ZG2 

Short  (S) 

4.0 

mm 

Medium  (M) 

8.5 

Long  (L) 

13.5 

16.5 

with  the  idea  that,  if  the  CCD's  with  high-speed  I/O  structures  were  too  last  to  determine 
an  upper  operating  speed  limit  for,  one  could  determine  the  limit  for  devices  with  I/O 
structures  designed  to  be  slow  and  extrapolate  to  the  performance  of  the  high-speed 
structures. 

The  fabrication  difficulties  described  above  took  their  toll  on  the  small  test 
devices  also.  Many  were  totally  defective;  many  operated  but  with  strange  characteris¬ 
tics;  and  all  were  highly  susceptible  to  gate  breakdown.  However,  some  operable  devices 
were  obtained  and  characterized.  After  the  data  was  analyzed,  a  surprisingly  close 
agreement  was  obtained  with  the  theoretical  results  presented  in  Section  2. 

5.2  Experimental  Results 

5.2.1  Procedures 

The  purpose  of  the  experiments  performed  on  these  two  devices  was  to 
investigate  the  filling  and  equilibration  (spilling)  phases  of  the  potential  equilibration 
input  method  as  a  function  of  the  time  allowed  for  each  phase.  This  was  done  while 
running  the  CCD's  at  a  modest  overall  clock  rate  of  1  MHz.  Thus  the  output  structures  of 
the  CCD  were  operating  well  within  their  performance  range  and  were  not  a  limiting 
factor. 

The  kind  of  clock  waveforms  used  are  shown  in  Figure  5-1.  The  width  and 
position  of  the  sampling  pulse  (SP)  applied  to  the  CCD  source  were  varied.  A  sample 
output  from  the  CCD  is  shown  in  Figure  5-2.  Feedthrough  of  the  various  operating  pulses 
can  be  seen  in  the  output.  During  the  first  100  ns  the  reset  pulse  (RP)  goes  high  and  back 
low,  O2  goes  high,  and  0j  goes  low,  causing  the  output  signal  to  appear.  To  each  side  of 
the  S00  ns  graticule  mark,  the  sampling  pulse  transitions  can  be  seen  (the  fill  time  is 
about  30  ns).  The  large  transition  just  past  900  ns  is  caused  by  capacitive  coupling 
between  0 2  and  the  floating  diffusion.  Thus  the  equilibration  time  is  thus  about  100  ns  in 
this  photo.  The  downward  ramp  in  the  signal  results  from  extremely  high  leakage, 
probably  in  the  reverse-biased  floating  diffusion. 
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5.2.2  Filling  Phase 


The  filling  phase,  as  expected  from  theory,  is  extremely  fast,  and  the  experi¬ 
ment.-..  results  are  not  easy  to  interpret  precisely.  No  effects  from  reducing  the  length  of 
the  sample  pulse  (SP)  could  be  seen  until  pulse  lengths  were  very  close  to  the  limit  set  by 
the  equipment. 

One  problem  was  to  determine  the  form  of  the  pulse  at  the  package  pin  of 
the  CCU.  The  only  available  tool,  an  FET  probe,  had  a  bandwidth  of  500  MHz.  Probe 
distortion  Is  snown  in  Figure  5-3.  Figure  5-3A  shows  SP  applied  to  the  50  ohm  input  of 
a  Tektronix  7104  oscilloscope  with  7A29  vertical  amplifier  plug-in  (BW  -  1  GHz).  Figure 
5-3B  snows  what  the  FET  probe  saw  when  SP  was  terminated  with  a  50  ohm  pad.  The 
apparent  pulse  amplitude  is  in  error.  The  apparent  pulse  length  is  correct,  however. 


i 
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Figure  5-3  -  Minimum  Length  Sample  Pulses 
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The  dependence  of  the  filling  behavior  on  SP  length  depends  critically 
on  the  exact  conditions.  Figure  5-4  shows  the  ouptut  signal  and  its  spectrum  for  a 
sine  wave  input  when  the  sample  pulse  caused  filling  to  occur  for  6  ns.  The  CCD  was 
of  type  L.  Almost  no  change  was  observed  as  the  fill  time  was  reduced  to  2.8  ns  (Figure 
5-5A).  At  2.4  ns  (Figure  5-5B)  distortion  can  be  seen.  The  oviput  is  unchanged  for 
small  signal  packets  (higher  output  voltage),  but  there  is  not  enough  time  to  provide 
more  than  a  half -full-well  when  the  input  signal  corresponds  to  a  large  pocket.  Some 
or  all  of  the  change  may  have  been  the  result  of  an  amplitude  change  in  SP  in  going 
from  2.8  ns  to  2.4  ns.  Reducing  the  SP  low  level  by  0.2  V  restored  the  signal  nearly 
to  its  original  form. 

Figures  5-6  and  5-7  show  another  set  of  results  with  pulse  lengths  of 
2.3  ns  arid  12  ns.  The  change  in  output  signal  is  so  small  that  it  can  be  seen  only  with 
the  spectrum  analyzer. 


Figure  5-4  -  Output  for  a  6  ns  Sampling  Pulse.  A)  Waveform;  B)  Spectrum. 
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Figure  5-7  -  Change  in  Result  when  SP  was  Increased  to  12  ns.  The  Sampling 
Pulse  is  Shown  in  A)  and  the  Spectrum  in  B) 

5.2.3  Equilibration  Phase 

The  dynamics  of  the  equilibration  phase  of  the  potential  equilibration  input 
technique  were  studied  by  using  a  fixed,  20  ns-long  sampling  pulse  which  ended  a  variable 
amount  of  time  before  0j  turned  on  and  transferred  the  charge  into  the  CCD.  Figure  5-8 
shows  results  for  a  medium  length  gate  CCD  with  equilibration  times  of  110  ns  and  5  ns. 
At  first  glance  the  pair  of  photos  seems  the  same  as  those  in  Figure  5-5.  Careful 
examination,  however,  reveals  a  critical  difference  between  the  two  distorted  outputs.  In 
Figure  5-8B  the  entire  output  signal  has  shifted  downward  (more  electrons  in  all  the 
charge  pockets)  until  saturation  of  the  CCD  transfer  wells  occurred.  In  Figure  5-5  the 
undistorted  part  of  the  output  is  unchanged.  The  explanation  is  that  in  Figure  5-8B  the 
input  well  under  IG2,  which  is  overfilled  during  the  sampling  pulse,  has  not  had  time  to 
equilibrate. 

At  this  point  a  very  important  result  was  observed.  A  slight  reduction  in 
the  DC  component  of  the  signal  on  IG2  caused  the  output  signal  to  be  restored  to 
its  original  form,  as  shown  in  Figure  5-9.  The  spectra  corresponding  to  the  signals  in 
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Figure  5-8A  and  5-9  are  shown  in  Figure  5-10;  they  are  virtually  identical.  The  most 
extreme  example  is  shown  in  Figure  5-11,  where  the  fill  time  was  2  ns  and  the 
equilibration  time  about  2  1/2  ns.  The  input  signal  bias  had  to  be  reduced  by  0.7  V  com¬ 
pared  to  the  value  for  long  fill  and  spill  times.  Thus  a  device  will  input  gates  8.5  U  m 
and  1 1.5  V  m  long  could  operate  using  a  potential  equilibration  input  at  a  frequency 
above  100  MHz.  Of  course,  full  equilibration  is  not  occurring,  and  some  of  the  advantages 
of  the  potential  equilibration  method  will  be  lost. 

In  order  to  get  a  quantitative  measure  of  the  speed  of  equilibration,  we 
measured  the  amount  of  input  signal  bias  shift  required  to  preserve  the  fidelity  of  the 
output  signal  as  the  equilibration  time  was  varied.  Plotting  the  results  on  semi-log  paper 
gave  remarkable  curves,  as  shown  in  Figure  5-12.  We  were  even  more  surprised  when  we 
replotted  the  theoretical  results  of  Chapter  2  on  semi-log  paper  and  found  the  same 
behavior. 
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Figure  3-10  -  Spectra  Corresponding  to  Figure  5-9B  (left)  and  5-10  (right) 
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Figure  5-1 1  -  Output  Signal  A)  and  Spectrum  B)  for  a  Medium  Gate  CCD 
with  a  2  ns  Long  Sample  Pulse  and  2  1/2  ns  of  Equilibration  Time 


Figure  5-12  -  Dependence  of  Input  Signal  Bias  Offset  on  Equilibration  Time 

5-10 


It  is  reasonable  to  suppose  that  the  slope  of  the  curve  depends  on  the 
amount  of  overfilling  during  the  fill  phase  of  the  input  cycle.  This  accounts  for  the  differ¬ 
ence  in  slopes  of  the  two  curves  in  Figure  5-12,  since  the  conditions  for  the  two  CCD's 
were  probably  not  the  same.  The  intercept  with  the  time  axis,  however,  should  be  a  good 
relative  measure  of  the  equilibration  time  for  the  structures.  The  theoretical  results  of 
chapter  2  indicated  that  the  equilibration  time  should  be  proportioned  to  the  length  of  the 
reference  gate  (10  j )  and  to  the  input  well  gate  (102)-  ^ie  values  of  these  lengths  for 
Figure  5-13  were  1pm  and  5  pm  respectively.  Applying  the  appropriate  factors  to  the 
1.6  ns  intercept  using  the  gate  lengths  from  Table  5-1  gives  intercepts  of  31  ns  and  71  ns 
for  the  medium  and  long  gate  CCD's  respectively,  very  close  to  the  experimental  values 
of  26  ns  and  69  ns. 


Figure  5-13  -  Excess  Input  Charge  Expressed  in  Volts 
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5.2.4  Conclusions 


Although  the  amount  of  data  obtained  was  limited  by  poor  device  yield 
resulting  from  processing  problems,  the  data,  not  only  was  in  qualitative  agreement  with 
the  theoretical  results  of  Chapter  2,  but  also  agreed  to  a  remarkably  quantitative  degree. 
There  is  little  doubt  that  the  ECMOS  input/output  devices,  had  they  been  fabricated  suc¬ 
cessfully,  would  have  been  capable  of  operation  at  frequencies  well  above  100  MHz. 
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